Nuclear constraints on non-Newtonian gravity at femtometer scale 
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Effects of the non-Newtonian gravity on properties of finite nuclei are studied by consistently 
incorporating both the direct and exchange contribution of the Yukawa potential in the Hartree- 
Fock approach using a well-tested Skyrme force for the strong interaction. It is shown for the first 
time that the strength of the Yukawa term in the non-Newtonian gravity is limited to log(|cv|) < 
1.75/[A(fm)] ' 5 +33.6 within the length scale of A = 1 — 10 fm in order for the calculated properties 
of finite nuclei not to be in conflict with accurate experimental data available. 
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The quest to unify gravity with other fundamental 
forces is among the most challenging scientific ques- 
tions for the new century [lj. While it is generally 
assumed that non-relativistic gravity obeys Newton's 
Inverse-Square-Law (ISL) for all distances greater than 
the Plank length of about 1.6 x 1CT 33 cm @-Q, the New- 
tonian gravity has been tested and various upper limits 
0-0] on the ISL violation have been set down to only 
about 10 fm so far Q. In fact, motivated by the possi- 
ble existence of both extra dimensions within string/M 
theories and new particles in the supersymmetric exten- 
sion of the Standard Model, it has long been proposed 
that the Newtonian gravitational potential between two 
objects of masses mi and mi at positions T\ and r-i may 
be modified to 191 
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where G is the gravitational constant, a and A are the 
strength and length scale of the non-Newtonian Yukawa 
potential, respectively Moreover, several modified grav- 
ity theories including the scalar-tensor- vector gravity [10( 
and f(R) gravity [1JJ lead to such a non-Newtonian po- 
tential at the weak-field limit. While the main pur- 
pose of this work is to constrain the parameters of the 
Yukawa term at femtometer scales, it is important to em- 
phasize that this kind of non-Newtonian potentials have 
been tested extensively against well-known observations 
at galaxy scales, see, e.g., Ref. [12[ for a recent review. In 
particular, we notice here that the non-Newtonian poten- 
tial of Eq. ([1]) has been used to successfully explain both 
the flattening of the galaxy rotation curves away from the 
Kepler limit, see, e.g. Refs. [lj, [lj] and the Bullet Clus- 
ter 1E0657-558 observations in the absence of dark mat- 
ter 1 151 1 . Moreover, it was shown that, unlike for massive 
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particles, the motion of massless particles is not affected 
by the Yukawa term. Consequently, all the lensing ob- 
servables obtained with the non-Newtonian potential of 
Eq. (|T|) are equal to the ones known from General Relativ- 
ity thanks to suitable cancelations in the post-Newtonian 
limit [lg,[l7|- Thus, to explain the gravitational lensing 
phenomenon, dark matter seems still needed with the 
non-Newtonian potential considered here. 

In the one-boson-exchangc picture, the Yukawa poten- 
tial may come from exchanging a light and weakly cou- 
pled spin-0 axion [18[ or spin-1 [/-boson [19[ correspond- 
ing respectively to an attractive or repulsive potential, 
and it can be further written as 



Vv(ri,r 2 ) = ± 
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where fi — 1/A is the boson mass, and g = -\/47r|ev|Gm 2 
is the boson-nucleon coupling constant with m being the 
nucleon mass. The exchanged boson may mediate the 
annihilation of dark matter particles [2fJ, [21J and it is 
related to several interesting new phenomena in particle 
physics and cosmology [22|, [23[ • It was also found recently 
that the Yukawa term affects significantly the Equation 
of State (EOS) of neutron-rich nuclear matter, thus prop- 
erties of neutron stars (NSs), e.g., their mass-radius rela- 
tion, moment of inertia and the core-crust transition den- 
sity/pressue [24| - |27| . While observed properties of NSs 
can be very well reproduced by combining the Yukawa 
term using appropriate parameters with nuclear EOSs 
that otherwise failed to do so by themselves [24|, [25( , the 
Yukawa term is not absolutely necessary because of the 
poorly known nature of the EOS of dense neutron-rich 
nuclear matter. On the other hand, it is not ruled out 
cither. To our best knowledge, no quantitative constraint 
exists on the strength of the Yukawa term at femtometer 
or shorter length scale. 

Nuclei are testing grounds of fundamental interac- 
tions among nucleons in the femtometer range. Gener- 
ally speaking, experimental measurements of global prop- 
erties of stable nuclei, such as their binding energies 
and charge radii, are very well reproduced consistently 



within about 2% uncertainty using established nuclear 
many-body theories and our current knowledge about 
the strong, weak and Coulomb forces without consider- 
ing gravity. Moreover, these properties are most strongly 
influenced by the well known isoscalar part of the nu- 
clear strong interaction around the saturation density of 
nuclear matter. Thus, in order to be consistent with 
current theoretical understanding and experimental ob- 
servations in nuclear physics, effects of gravity on prop- 
erties of stable nuclei have to be less than about 2%. Us- 
ing this requirement and by consistently incorporating 
both the direct and exchange contributions of gravity in 
the Hartree-Fock (HF) approach using a well established 
Skyrme force for the strong interaction, we show for the 
first time that the strength of the Yukawa term is lim- 
ited to log(H) < 1.75/[A(fm)] - 54 + 33.6 in the range of 
A = 1 — 10 fm. This provides a reliable reference to test 
stringently non-Newtonian gravitational theories in this 
previously unexplored region. 

The contribution from the non-Newtonian Yukawa po- 
tential to the potential energy of the nuclear system is 
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where P r , P a , and P T are the space, spin, and isospin 
exchange operator, respectively, and \i) is the quantum 
state of the «th particle containing spatial, spin, and 
isospin parts. The first term in Eq. ([3]) is the direct con- 
tribution and can be calculated from 
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where p(r) = J2i. a T ^Tii 1 ?, (r)</>Ti(r, a) is the nucleon 
density with cj) T i(r,cr) being the spatial wave function 
of the ith particle with spin a and isospin r. Using 
P a = (1 + (?i • a 2 )/2, with 0V2) being the Pauli oper- 
ator acting on the first (second) term, the second term 
in Eq. © representing the exchange contribution can be 
expressed as 
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where p T (ri,r 2 ) = YLi,a < t > %i^ f i^)4>Ti{r-z,cr) and 
Pr(ri,r 2 ) = T lt T l a,a' ( l } U( r 'i^')4'T t (r 2 ,cr){a'\a\a) 
are the off-diagonal scalar and vector part of the 
density-matrix, respectively. Introducing the coordinate 
transformation r = [r\ + r 2 )/2 and s — f\ — r 2 , they 
can be calculated from the density-matrix expansion 
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In the above, k T = (3ir 2 p T ) 1 is the Fermi momentum, 
r T (r) = y*. |V</) T i(f, a)\ 2 is the kinetic energy density, 

and J T {r) = -iEiE CT ,<7' Kii?, <r')V<l>ri{r,cr) x (a'\a\a) 
is the spin density. The psL(k T s) and g(k T s) can be ex- 
pressed in terms of the first- and the third-order spher- 
ical Bcsscl function as psL{k T s) — 3ji(k T s)/(k T s) and 
g{k T s) = 35j3(/e T s)/[2(fc T s) 3 ], respectively. 

Using the expressions from the density-matrix expan- 
sion, the potential energy density from the exchange con- 
tribution of the Yukawa potential can thus be written as 

B$(f) = ^{A[p T (r)] + J B[p T (r)]r r (r) 

t— n,p 

+ C[ Pr (f)][Vp r (f)] 2 + ^[p r (f)]J r 2 (r)}, (8) 
where 
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with V Y {s) =g 2 e-^ s /(4:TTs). 

To calculate properties of finite nuclei within the HF 
approach, we solve the Schrodingcr equation 






(13) 



where m*, U T and W T are the nucleon effective mass, the 
single-particle potential, and the form factor of the one- 
body spin-orbit potential, respectively. From the varia- 
tional principle, the Yukawa contribution to the single- 
particle potential has a direct term U YU and an exchange 



term U^ E , namely, 
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The Yukawa contribution to the form factor Wj is 

W? = 2<p(p T )J T . (17) 

In addition, the effective mass is modified by the Yukawa 
term according to 
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Most available nuclear effective interactions have been 
tuned to reproduce not only empirical properties of sym- 
metric nuclear matter at saturation density but also 
global properties of finite nuclei along the j3 stability 
line. These interactions differ mostly in their isovector 
parts and/or the effective three-body forces used. In 
this work, we use properties, such as the charge radii 
and binding energies of stable medium to heavy nuclei 
that are not much affected by the still uncertain parts of 
the strong nuclear interaction. Specifically, we use here 
the MSLO parameter set of a Skyrme-like force which 
has been tested extensively and was shown to describe 
the charge radii and binding energies of finite nuclei very 
well [Ij. 

Before exploring effects of gravity on properties of fi- 
nite nuclei, it is instructive to first compare the magni- 
tudes of the Coulomb and Yukawa potentials between two 
protons. Shown in Fig.[TJare the Yukawa potentials using 
a = -1.24 x 10 36 or -1.24 x 10 34 , and A = 1 fm or 10 fm, 
respectively. With a = —1.24 x 10 36 , which leads to the 
coupling constant c/ 2 /47r = 1/137, and a long interaction 
range of A = 10 fm, the Yukawa potential approaches the 
Coulomb potential at short distances, while it decreases 
faster towards long distances for smaller values of A. 

As an example illustrating effects of the Yukawa po- 
tential on properties of finite nuclei, shown in Fig. [2] are 
the charge density profiles of 208 Pb obtained within the 
HF approach using the MSLO interaction with different 
Yukawa parameters. With a on the order of 10 34 , the 
Yukawa potential has negligible effects on the charge den- 
sity profile. Increasing the magnitude of a to 10 36 , effects 
of the Yukawa potential become significant especially for 
larger values of A. 

Next we investigate more systematically how the 
Yukawa potential affects the charge radii r c and bind- 
ing energies (B.E.) of three typical medium to heavy 
nuclei by varying continuously the strength parameter 
a. Shown in Fig. [3] arc the results for 208 Pb, 120 Sn, 



10 1 



> 10° 

0) 



> 

10" 2 
10" 3 



1 ' i ' i 

Coulomb potential 

— - a = -1.24x10 36 , X = 10 fm 

a = -1.24x10 36 , X= 1 fm 

a = -1.24x10 34 , X = 10fm 




-+ 



r(fm) 



FIG. 1: (color online) A comparison of the Coulomb and 
Yukawa potentials between two protons as a function of sep- 
aration with different parameters. 
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FIG. 2: (color online) Charge density profiles of 208 Pb from 
MSLO with and without the Yukawa potential from different 
parameters. 



and 48 Ca with A = 1 fm and 10 fm, respectively. The 
horizontal lines are the mean values of the experimental 
data [3l|, [32| ■ As the Yukawa potential varies from be- 
ing repulsive (with negative values for a) to attractive, 
the charge radius decreases and the nuclei becomes more 
bound. The effect is roughly proportional to the value of 
a. Moreover, the effects arc larger for heavier nuclei due 
to the finite-range nature of the Yukawa potential. 

As mentioned earlier, essentially all existing nuclear ef- 
fective interactions can describe the experimental charge 
radii and binding energies of medium to heavy nuclei 
within about 2%. The latter is thus the largest room 
available to accommodate effects of gravity. As shown 
by the green bands in Fig. [3] this then allows us to set 
a nuclear upper limit on the strength of the Yukawa po- 
tential in the length range of 1 to 10 fm. Here we use the 
Yukawa potential effects on 208 Pb which is more strongly 




40 -8 -4 

ot(10 35 ) 



FIG. 3: (color online) Charge radius r c [(a), (b)] and binding 
energy (B.E.) [(c), (d)] of 208 Pb, 120 Sn, and 48 Ca as functions 
of a at different length scales. The results at a = are 
those from the MSLO interaction only. The green bands are 
constrained ct values from taking 2% as the uncertainty of the 
charge radii and binding energies. 



affected by gravity. More specifically, the constraint from 
the binding energies can be parameterized as log(|a|) < 
1.75/[A(fm)]°- 54 + 33.6 while that from the charge radii 
can be written as log(|a|) < 1.18/[A(fm)] 079 + 35.0 for 
A = 1 — 10 fm. Since the binding energy is more sensi- 
tive to the strength of Yukawa potential, the first con- 
straint is more stringent and can be used as the nu- 
clear upper limit. This constraint thus limits the cou- 
pling constant of the exchanged bosons with nucleons 
to log( ff 2 ) < 0.10[/i(MeV)] 054 - 3.53 if the mass of the 
bosons is between 20 MeV and 200 MeV. Previously, the 
shortest range probed is above 10 fm in neutron scatter- 
ing experiments [8j, [33j-|36j . For a comparison, shown in 
Fig. 2] arc the various upper limits on the magnitude of 
a in the range of 10~ 15 to 10 -10 m. The constraint ex- 
tracted from the present work extends to the previously 
unexplored region of 1 — 10 fm. Overall, there is a clear 
trend that the allowed deviation from Newtonian grav- 
ity increases as the interaction range decreases, reflecting 
the increasing difficulties of the measurement. 

In summary, effects of the non-Newtonian grav- 
ity on properties of finite nuclei were studied within 
the Hartree-Fock approach incorporating a well-tested 
Skyrme force for the strong interaction and the non- 
Newtonian gravitational potential. For the first time, 
the strength of the Yukawa term is limited to log(|a|) < 



1.75/[A(fm)] a54 + 33.6 within the length range of A = 
1 — 10 fm in order for the calculated properties of finite 
nuclei not to be in conflict with the very accurate exper- 
imental data available. This constraint serves as a use- 
ful reference in constraining properties of weakly-coupled 
gauge bosons and further explorations of possible extra 
dimensions at femtomcter scale. 
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(Barbieri et a/. 1975} 

(b) n-mirror scattering 
(Nesvizhevsky et al. 2004) 
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FIG. 4: (color online) The constraints of a for A = 1 — 10 
fm from charge radii and binding energies of heavy nuclei and 
those at longer distances extracted from analyzing neutron 
scattering experiments in Refs. [a. |33B36| . 
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